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Abstract 



We study the long-time behavior of an exterior Hele-Shaw problem in random media 
with a free boundary velocity that depends on position in dimensions n > 2. A natural 
rescaling of solutions that is compatible with the evolution of the free boundary leads 
to homogenization of the free boundary velocity. By studying a limit obstacle problem 
for a Hele-Shaw system with a point source, we are able to show uniform convergence of 
the rescaled solution to a self-similar limit profile and we deduce that the rescaled free 
boundary uniformly approaches a sphere. 



We consider a Hele-Shaw type problem (HS) in random media and study long-time 
behavior of the solutions. Let n > 2 and let K C Oo C M. n , K be a nonempty compact 
set and f^o be a bounded open set such that K and Clo have smooth boundaries (see 
Figure HJ. Let T = dK and Tq = 8{Iq. Formally, the Hele-Shaw- type problem is to find 
a function v(x,t) : M. n x [0, oo) — > R that solves 
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Figure 1: Initial configuration 

Here Dv and Av are, respectively, the gradient and Laplacian of v, with respect to the 
space variable x, and vt is the partial derivative of v with respect to the time variable t. 

Let (A, F, fj) be a probability space. We consider g(x,u) : W 1 xvl-^Ka continuous 
function in x for a.e. uj £ A, satisfying 

< m < g{x,u) < M for all x 6 M n , a.e. ui € A, (1.1) 

where m and M are positive constants. Furthermore, in order to observe some averaging 
behavior as t — >■ oo, we assume that g is stationary ergodic. In other words, we assume 
that for every x E M n there exists an ergodic transformation t x : A^t A such that 

g(x + x' , uS) = g(x' , t x uj) for all x, x € M n and a.e. w G A. (1-2) 

We say that a measure preserving transformation r x is ergodic if the only invariant sets 
are either zero or full measure, that is, if B C A such that t x (B) = B for all x £ M n , 
then [i{B) = or 1. For a more detailed discussion on stationary ergodic media, see for 
instance [HITT]. 

The classical Hele-Shaw problem in two dimensions with g = 1 was introduced in |6J 
modelling a slow movement of a viscous fluid injected in between two parallel plates close 
to each other that form the so called Hele-Shaw cell. This problem naturally generalizes 
to all dimensions n > 1. 

We study a Hele-Shaw type problem (|HS|) describing a pressure-driven motion of 
a fluid in an inhomogeneous random medium where the velocity law of the fluid on 
the free boundary depends on the position. The homogenization of this problem was 
recently studied in [8]. Free boundary problems with similar velocity laws have various 
applications in the plastics industry p3], in electromechanical machining [10] and others. 
In fact, Hele-Shaw problem can be thought of as a quasi-stationary limit of the Stefan 
problem with a similar boundary velocity law, modelling heat transfer, see |12yi6j. 

In this paper, we are concerned with the behavior of the solution for large times. 
The main result, Theorem EU is the locally uniform convergence of a rescaled viscosity 
solution of the inhomogeneous problem (jHS]), 

v x {x, t) = A (n - 2)/n w(A 1/n x, Xt), if n > 3, 

(see section [3] for discussion of the rescaling and the case n = 2), which then satisfies 



v £ = g (\ 1 / n x,uj) Dv x 2 on d<|V > o} 
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and its free boundary as A — > oo to the self-similar solution of a Hele-Shaw problem 
with a point source, formally 

-Av = C5 in {v > 0}, 

1 



IV ■ ' ' (1-3) 



vt = -j—^-\Dv\ on d{v > 0}, 

V(;0) =0 

Here 5 is the Dirac J-function with mass at the origin and the constant C depends on 
K and n. Constant /|y represents the "average" of | and it will be properly defined 
later. 

The following heuristics can help with understanding the boundary velocity of the 
limit problem. For the moment, we suppose that g is periodic and the solution v is 
sufficiently smooth. Since the free boundary d{v > 0} is a level set of v, its normal 
velocity is given by V = r^b- If we interpret (|HSp as a quasi-stationary limit of the 
one-phase Stefan problem where v is the temperature and —Dv is the heat flux, the 
quantity ~ represents the latent heat of phase transition depending on position (see 
[12]). Indeed, for the free boundary to expand by dx, the heat flux must deliver -dx of 
heat into the interface. Since it delivers exactly \Dv \ dt of heat in time dt, we have 

dx i 

V~ — = —\Dv\=g(x)\Dv\. (1.4) 

When we rescale the problem, the latent heat will simply average out and therefore 

V = -^-\Dv\ 

is the natural free boundary velocity of the limit solution. 

A similar result to ours was obtained in [13] for weak solutions of the homogeneous 
Hele-Shaw problem with free boundary velocity 

V = \Dv\ . 

In the current situation, however, the velocity law of the free boundary depends on 
the position and therefore the techniques from |13] can provide only lower and upper 
bounds on the free boundary radius. This requires us to use a more refined method to 
prove the convergence of the solution to the self-similar asymptotic profile. We combine 
the strengths of two notions of solutions of (jHS]) - viscosity and weak - using their 
correspondence studied in [Sj. 

The weak solution u(x, t) : W 1 x [0, oo) — > R was introduced in [3] by applying the 
Baiocchi transform 

u(x, t) = v(x, s) ds. 
Jo 

The function u formally solves the Euler-Lagrange equation 

-An = — 7-^Xm.™\n in {u > 0}, 

u = \Du\ = ond{u>0}, ^' 5 ^ 

u = t on K 
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of some obstacle problem. Since g satisfies (jl.ip and (|1.2p . the subadditive ergodic the- 



as A — > oo. This then leads to homogenization of the obstacle problem (|1.5p . 

While in [8] the uniform convergence of the weak solutions in the homogenization 
limit follows more directly from the regularity of the obstacle problem (j 1 . 5 j) , the current 
situation is considerably different. As the fixed domain K shrinks to the origin due to 
the rescaling, the rescaled weak solution gets singular. Thus the convergence of weak 
solutions as A —7- oo is restricted to a locally uniform convergence on compact sets not 
containing the origin and the obstacle problem regularity cannot be used directly. We 
resolve this difficulty by introducing a limit obstacle problem and studying its wellposed- 
ness. The main challenge is posed by the singular behavior of the solution as \x\ — > 0. 
From this development, we can deduce the necessary local uniform convergence away 
from the origin. As the last step, we extend the applicability of the tools developed in 
[8] to our situation with a singularity at the origin. This allows us to conclude that 
the rescaled free boundary converges uniformly to the free boundary of the asymptotic 
profile, a sphere. 

The current method requires g(x, uS) to be strictly positive, the main reason being the 
loss of uniqueness of viscosity solutions if g becomes zero and the free boundary can stop. 
It is an interesting open problem if g can reach zero yet | 6 L\ oc so that (II. 6p still makes 
sense. It is currently under investigation by the author. The interpretation (jl.4p of (|HSp 
in terms of heat balance suggests that the homogenized velocity law V = (l/g}~ \Dv\ 
continues to hold, at least for maximal viscosity solutions. 

A brief outline follows: Section 2 summarizes the two notions of solutions for problem 
(IHS|) and their correspondence. In section 3, we introduce the rescaling of (|HSp . which 
is then followed by an investigation of radially symmetric test functions of section 4. 
The main contribution of the paper is then contained in sections 5 and 6, where we 
first define and study the limit obstacle problem and then show that the rescaled weak 
solutions of (IHSp converge to the unique solution of the limit obstacle problem. Finally, 
equipped with those results in section 7, we turn back to the viscosity notion to prove 
the uniform convergence of free boundaries to those of the limit problem. 



Notation. Since the arguments are independent of to, we fix to E A for which all ([Lip . 
(II. 2p and (II. 6p hold and we omit u in the rest of the paper. 

Throughout the article, we will make use of the standard notation for the bilinear 
form •) on H 1 ^) and the scalar product (•, -}q on L 2 (Q) for some domain J7, 





for a.e. uj £ A 



(1.6) 



2 Hele-Shaw type problem 




Du ■ Dv dx 




uv dx. 



Whenever = M. n , we drop the subscript f2. 
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For a set A, we denote A c its complement. Given a nonnegative function v, we 
respectively define the positivity set of v and the free boundary of v, 

fi(u) := {(x,t) : v(x,t) > 0}, T(v) := dn(v), 

and their time-slices 

fit(u) := {x : u(ar, i) > 0}, T t {v) := dQ t (v). 

(/)+ will be the positive part of /, i.e. 

(/)+ = max{/,0}. 

Lastly, B r (x) is the open space ball of radius r centered at x, 

B r (x) = {y € R n : \y - x\ < r}. 

This section reviews the two notion of solutions of (jHSj) that will be used in this 
paper. 

The notion of viscosity solutions for (|HS|) was introduced in [7] and used in [8] to 
prove uniform convergence of the free boundary in the homogenization limit. We will 
consider solutions on the spacetime cylinder Q = (R n \ K) x [0, oo). 

Definition 2.1. A nonnegative upper-semicontinuous function v defined in Q is a vis- 
cosity subsolution of (|HS|) if 

(a) for each T > 0, the set £l(v ) n {t < T} n Q is bounded, and 

(b) for every cj) S C^.'\{Q) such that v — <fi has a local maximum in Q(v) n {t < to} H Q 
at (xo,to), the following two conditions hold: 

(i) if v(xo,to) > then — A(j)(xo,to) < 0, 

(ii) if (aJo,*o) ^ |jD0| (xo,to) / and — A0(xo,to) > 0, then 

(ci>t-g(xa)\D<f>\ 2 )(xo,ta) < 0. 

Definition 2.2. A nonnegative lower-semicontinuous function v defined in Q is a vis- 
cosity supersolution of (jHSj) if for every <fi £ C^' t (Q) such that u — has a local minimum 
in n {t < to} H Q at (xo,to)j the following two conditions hold: 

(a) if v(xo,to) > then — A0(xo,to) > 0, 

(b) if (x ,to) S r(«), |D0| (x ,t ) / and -A</>(x ,t ) < 0, then 

(0 i - 5 (x o )|^| 2 )(xo,to) > 0. 

Now we can define viscosity sub- and supersolutions with appropriate initial and 
boundary data. 

Definition 2.3. A viscosity subsolution v of (IHSp in Q is a viscosity subsolution of 
(|HS|) in Q wt/i initial data vo and boundary data 1 if 

(a) v is upper-semicontinuous in Q, v = vq at t = and u < 1 on T, 

(b) Q(v) n {i = 0} = {x : u (x) > 0}. 



6 



N. Pozar 



Definition 2.4. A viscosity supersolution v of (jHS]) in Q is a viscosity supersolution of 
(|HS|) in Q with initial data vq and boundary data 1 if 

(a) v is lower-semicontinuous in Q, v = vq at t = and u > 1 on T, 

Finally, we can define a viscosity solution: 

Definition 2.5. A viscosity supersolution v of (IHSp in Q (with initial data vq and 
boundary data 1) is a viscosity solution of (|HS|) in Q (with initial data vq and boundary 
data 1) if 

v*(x, t) := limsup v(y,s) 
(y,a)-+(x,t) 

is a viscosity subsolution of (|HS|) in Q (with initial data vq and boundary data 1). 

The following comparison principle is due to [7J: 

Theorem 2.6 ([7J Theorem 1.7]). Let x>\ and v% be respectively viscosity sub- and 
supersolution of (|HS|) . If {vi(-,0) > 0} is bounded, v\{x, 0) < i>2(x,0) in {ui(-,0) > 0} 
and v\ < V2 on dK x [0, oo), then 

Vi < v 2 on K c x [0, oo). 

The second notion of solutions is the one introduced by [3], given by an obstacle 
problem for each time for a new variable u(x, t) = f Q v(x, s) ds: 

Definition 2.7. The function u(x, t) is called the weak solution of problem (|HSp if for 
every t > 0, w = u(-,t) solves the obstacle problem 

w e JC t , 

( 2 1 ) 

a(w,ip - w) > (--^jXR"\n ,¥ ~ w) for all cp E Kf, 
where 

Kt = {<p G ^d( Mn ) 5 <P > in n > <^ = t on K}. 

The theory of the obstacle problem (|2.ip is well understood. We refer the reader to 
|15j for instance. We have the following: 

Proposition 2.8. Suppose that K C Oo C IR" - , Oo is bounded open set and K is a 
compact set with G C > . Then f)2. 1 1) has a unique solution u(-,t) G ^(M™) for all 
t > 0. For every T > there exists -R(T) such that 

n t (u) C Sfl(0) for all t G [0, T]. 

Moreover, for any 1 < p < oo, there is a constant C = C(m,n,p) such that 

[K- ) *)llw r 2.p(R»\iC) - ^ for a11 * - °- 
Finally, the function 1 1— )• u(x,t) is Lipschitz continuous with 

|u(x, ti) — u(x, *2 ) | < |*i — t 2 \ for all < t\ < t 2 , 

for all x6l n 
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One of the main tools used in this paper is the following correspondence of weak and 
viscosity solutions proved in [8]: 

Theorem 2.9 (cf. [8j Theorem 3.1]). Let u(x,t) be the unique solution of (12. ip . For 

each t > 0, let v(-,t) be the solution of 

'Av(-,t) = in fit (it) \ K, 
< v = 1 on K 

v = T t (u), 

i.e. for every t > define v(-,t) as the supremum of all lower semicontinuous functions 
w(x) for which there exists some s, < s < t, such that 

— Aw < in O s (it), 

< w = 1 on if, 

supp w C O s (xt). 

Then v is a viscosity solution of (|HS|h Moreover, 

t) = d^~u(x, t), 

the left time derivative of it. 



3 Rescaling 

We are interested in the behavior of (IHSj) as t — > oo. The goal is to prove that a 
rescaled solution converges to the self-similar solution of the homogeneous Hele-Shaw 
problem with a point delta source. As was observed in [13], the behavior depends on the 
dimension and the necessary rescaling is more involved in the two-dimensional case. 



3.1 Case n > 3 



The heat interpretation f j 1 . 4 [) of the free boundary velocity from the introduction sug- 
gests that the radius of the free boundary should behave as ~ t 1 ' 71 and thus we use the 
rescaling 



If we define 
(IHSD : 



v x (x, t) = A( n " 2 )/ n i;(A 1 / n x, Xt). (3.1) 
-[7^- and Oq := ^r^r, we derive the following problem for v x from 



-Ai/ 







; \ = A (n-2)/r. 



in n(v x )\K x , 
on K x , 



v x = g x (x)\Dv x \ on T(v x ), 
v x (-,0)=v x onn x \K x , 

where g x {x) = g(X 1 ^ n x) as can be easily seen from 

X (n-2yn+l Vt ( X l/n x ^ ^ = g X ^ A 2(n-2)/n+2/n ^ Dv ^\ l / n X, Xt) 

v t (X 1/n x,Xt) = g x (x) (Dv)(X 1/n x,Xt) 



(3.2) 



(3.3) 
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compared to the boundary condition in (|HS|) 

v t {X 1/n x,Xt) = g(X 1/n x) (Dv){X l ' n x,Xt) 

The rescaled problem (13.2p implies the weak formulation for the rescaled u x (-,t) 
JqV X (-,s) ds, the obstacle problem 



a(w,cp- w) > (-jxXmrv\ n x,<p- wj for all if G /C A , 



(3.4) 



for each t > 0, where 



K, x = i^p G ff^R"), > 0, = A (n ~ 2)/n t on K x }. 



The relation of u x to u can be deduced from the definition of u x as 

u x (x,t) = [ v x (x,s) ds 
Jo 



A (n-2)/n t l v ^X l l n x,\ S ) ds 

Jo (3.5) 



v da 



= A (n-2)/n / u(A l/n X)(T) ^ 

Jo X 

= A- 2 / n n(A 1 / n x,At). 

The following observation will be useful: 

Lemma 3.1. The harmonic function p(x) = C\x\ 2 ~ n is invariant under the rescaling 
(|3.ip . In particular, if v(x,t) < C|x| 2_n for all i > then also 



v x (x,t) < C\x\ 2 ~ n foralli>0. 

Proof. Observe that 

p x ( x ) = CX^ n - 2 ^ n \X^ n x\ 2 - n = C\x\ 2 - n . 

□ 

3.2 Case n = 2 

Since the rescaling (|3.ip does not recover the logarithmic singularity at the origin in 
dimension n = 2, we use a different rescaling that preserves the system: 

v x (x,t) = logK(X)v(K(X)x,Xt), 

where TZ(X) is the unique solution of 

7l 2 log ft = A (3.6) 
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for any A > 0. In fact, the solution can be found explicitly in terms of the Lambert 
W function, 1Z(\) = exp (iW(2A)). It is easy to see that lim,x_ ) . 00 7£(A) = oo. Taking a 
logarithm of (|3.6p gives 



21ogft(A) 1 + 



log log log A 



We deduce 

which together with f|3.6[) yields 



21ogft(A) 21ogft(A) 
log A 



0. 



lim 

A^oo 2 log 7£(A) 



lim 



A— >oo 7^.oo(A) 

In this case, we define if A : 



^oo(A) 



2A 
log A 



1/2 



(3.7) 



(3.8) 



^y. Since v satisfied (IHSp . u A satisfies 



-Au A = 



logft(A) 
g A (x) \Dv 



Al 2 



V A (-,0) 



" W(A) 

in ft(u A ) \ if A , 
on K A , 
on T(v A ), 
on ft A \ if A . 



(3.9) 



Here g x (x) = g(lZ(\)x), which can be derived as in (13.3H together with (|3.6j) . 

The rescaled problem (|3.9j) implies a weak formulation for the rescaled u x (x,t) 
Jq v x (x, s) ds of the form (|3.4p where now 

/C A = |(^ E il^M"), 9? > 0, = logft(A)t on if A }. 

This induces scaling for weak solution u analogous to (|3.5p . 



u A (M) = l ° g ^ X \ (K{\)x,\t). 
A 

Remark 3.2. Note that the scaling in dimensions n > 3 and the scaling in dimension 
n = 2 are qualitatively different. Indeed, in dimension n > 3 a solution with a point 
source and a solution with a finite source have the same asymptotic speed of the free 
boundary ~ t 1 / 11 . That is not true in dimension n = 2 where a point source solution 

has asymptotic speed ~ t 1//2 , but a solution with a finite source has a slower asymptotic 

/ \V 2 
s P eed ~ (&t 



4 Comparison with radially symmetric solutions 

Radially symmetric solutions for the Hele-Shaw problem, derived in [13J, will serve as 
test functions in our arguments. A radially symmetric solution in the domain |x| > a, 
t > is a pair of functions p(x, t) and R(t), where p is of the form 

Aa 2 ~ n (|x| 2 ~ n - R 2 ~ n (t)) 
v(x,t) = s if n > 3, (4.1a) 
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or 

p(x,t) = — J^T^ ifn = 2, (4.1b) 

and .R(i) satisfies a certain algebraic equation (see [13] for details, we will be interested 
only in the behavior as t — > oo). This solution satisfies the boundary conditions 



(4.2) 



p(x,t) 


= Aa 2 ~ n 


for |x| 


= a > 0, 


p(x,t) 


= 


for \x\ 


= R(t), 


R'(t) 


= \\D P \, 


for \x\ 


= R(t), 


R(0) 


= b> a. 







Also 



or 



R(t) {An(n-2)\ 1/n c tt . 

lim — ^- = 1, Coo = ^= >- if n > 3, (4.3a) 

Coot 1 /" \ L I 



lim ^ 1/9 =1, Coo = 2 V / I7^ if n = 2. (4.3b) 



In dimension n = 2, we will also make use of the limit behavior 

lim — — — = -. (4.4) 
t^oo logt 2 

Lemma 4.1. For L = ~ (resp. L = -A*), with m, M defined in (11. lj) . i) is a viscosity 
subsolution (resp. supersolution) of (|HS|) in Q. 

Proof. Assume that = m. The case = M is similar. We will show that p(x, t) 
satisfied the conditions in definition 12.11 

p(x,t) is clearly continuous and (a) immediately follows from (I4.1a|) and (14.31) . 

As for (b), (i) follows from the fact that p(-,t) is harmonic in f?t(p) for every t > 0. 

To show (ii), let <f> be a smooth function such that p — 4> has a local maximum in 
fi(p) n {t < t } at (xo,*o) G r(p) such that |Z)0| (a; , *o) / and -A0(x o ,to) > 0. Let 
-B r C ^ (p) be a ball touching I\ (p) from inside at (xo,to) in which —Acfr > (by 
smoothness of boudaries and continuity of A0 there is such a ball). Then by the Hopf's 
Lemma (see [4], ch. 6.4.2), since — A(</> — p) = —A(p < 0, 

P) > at(x ,t ), 

ov 

where v is the inner normal. Thus together with the fact that the boundaries are smooth 
level sets, we conclude that 

dp dip 

\Dp\(x ,t ) = 7^ < = \D<p\(x ,t ). 
The normal velocity of the free boundary is given as V n = pt/\Dp\ (level sets), we have 

lH r(aJo,*o) > tttttO^o^o)- 
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Figure 2: Arrangement of domains for comparison 



Hence finally 



(^-<?(* o )|Z)0| 2 )(xo,to) < 1^1 \ m^-9( x o)\ D P\) (*o,to) 

®<\D<j>\ (^- -m\Dp\^ (x ,t )^0. 



□ 



Comparison with the radially symmetric sub- and supersolution provides us with the 
following straightforward result: 



Lemma 4.2. Let v be a viscosity solution of (jHSJ). There exists to > and constants 
pi and p2, < pi < p2, such that 

P\t l l n < min Id < max Id < pot 1 / 71 if n > 3, 

r t (v) r t (v) r 

or 

piR.it) < min Id < max Id < p2R.it) if n = 2, 

r t (t>) r t (u) 



for t > to and 



max x < P2 

r«(«) 



for < t < t . 

Moreover there is a constant C > such that 

< v(x,t) < C\x 



2-n 



Proof. Given problem (|HS|) . we can find constants a, b± such that B a (0) C K and 
B b _ (0) C ft C C B b+ (0) (see Fig. EJ). Set also L_ = L+ = jj. Using the 
maximum principle for harmonic functions, we can also find constants A- and A + so 
that p- and p + , the radially symmetric subsolution and supersolution, respectively, 
given in (|4.ip that satisfy the boundary conditions (|4.2|) with the respective constants, 
also satisfy 

p_ -<v ^. p + on (ivT c x {0}) U (r x (0, oo)) . 
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Then by the comparison principle for viscosity solutions, Theorem 2.6 in [5J, we have 

P- -< v -< p + on K c x (0,oo). 

and also 

v(x,t) < p+(x,t) < C\x\ 2 ~ n , 

for C > large enough and 

R-(t) < min \x\ < maxlxl < R+(t), 
r t (v) r t (v) 

where R± are the radii of free boundaries of functions p±, respectively, and their limit 
behavior is given by (|4.3p . □ 

Lemma 14.21 is all we need to prove the crucial theorem from |13j : 

Theorem 4.3 (Near-field limit). The viscosity solution v(x,t) of the Hele-Shaw-like 
problem (|HS|) converges to the unique solution P{x) of the exterior Dirichlet problem 

AP = 0, x G R" \ i^, 

P = 1 - * 61 "- (4.5) 

lim |x|->oo p ( x ) = °> if « > 3, or 

P is bounded if n = 2. 

as £ — )• 00 uniformly on compact subsets of K c . 

Proof. See proof of |13^ Theorem 4.1]. □ 
The following constant C* characterizes the singularity of the limit solution: 

Lemma 4.4. There exists a constant C* = C*(K) such that the solution P of problem 
(14. 5p satisfies 

lim |x| n - 2 P(x) = C*. 

|cc|— >oo 

Proof. See Lemma 4.5 in [13J. □ 

Finally, we also need to improve on the convergence result of radially symmetric 
solutions from 1131: 



Lemma 4.5. Let p(x, t) be a radially symmetric solution (|4.ip of the Hele-Shaw problem 
satisfying the boundary conditions (|4.2p with constants A, a, b and L. Then the rescaled 
solutions p x (x, t) converge locally uniformly on the set (M. n \ {0}) x [0, 00) to the solution 
of the Hele-Shaw problem with a point source, 



V(x,t) = V A , L (x,t) = 

where 

Pit) = PL® 



A [\x\ 2 ~ n - p 2 - n (t)J , n>3 
A(log^) n = 2, 




An(n-2)t 
L 



1/n 

, n > 3 
n = 2 
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Proof. We shall show uniform convergence on sets {(x, t) : \x\ > e, < t < T} for some 
e, T > 0. Let t = p _1 (e/2). Starting with n > 3, we have 



Aa 2 - (|x| 2 ^ - (^t 1 /*) 
p X (x,t) = 



2-n 



a 2-n _ R2-n(Xt) 

There are two cases: 
(a) t <t<T 

Due to (|4.3j) . we have 

i?(At) _ i?(At)_ fl/n ^ / Ao"- 2 n(n-2)A 1/n = ^ 



A 1 /" (At) 1 /" ^ L 

R 2 ~ n (\t) -> 0, 

as A — >• oo uniformly on to < i < ^- That shows uniform convergence p x — > V on 
{M > e, t <t< T}. 

(b) < t < t 

Clearly V(x,t) = in {|x| > e, < t < t }. Since 

R(Xt) R(Xt ) ^,..,e 

for all A large enough, we see that p x = = V on {\x\ > e, < t < to} for all A 
large. 

In dimension n = 2, the rescaling yields 

(log 



Tl(X)\x\ 

P A (x,t)=A 



log - 



log^(A) 

Again, we split the proof in two cases: 
(a) t <t<T 

Rewriting R(\t)/lZ(\) and using (|4.3b|) and (|3.8j) . we obtain 

xt V/ 2 ( x W 2 



R(Xt) _ R(Xt) [T^xt) [TEix) 2 /A \ 1/2 tl/2 l_ 



yiogxt) \\ogxj 

as A — 7- oo uniformly for to < t <T. 
Similarly, flHJ) and fl37F|) lead to 

log^M _ bg jZ(At)- logo log At logA 1 
logft(A) log At ' log A ' logft(A) 2' 

as A — )• 00 uniformly for to <t < T. This proves uniform convergence of p A — >• V 
on {|x| > e, t < t < T}. 
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(b) < t < t 

Argue as in the case n > 3, (b). 

□ 



5 The limit problem 



Our current task is a characterization of the limit of rescaled weak solutions u x as 
A — > oo. We want to show that the limit satisfies a certain obstacle problem that can 
be interpreted as a Hele-Shaw problem with a point source. Existence and uniqueness 
of such a problem in 2 dimensions was studied in [2J. Our situation requires extending 
the definition to all dimensions n > 3. 

First define U Aj L{x,t) to be the Baiocchi transform of V Ai L(x,t), introduced in 
Lemma 14.51 (see proof of Theorem 15.11 for derivation) : 



U A ,L(x,t) = 

c \2~n + L_ N 2 _ 1 {Ant) 2/n („-2)(2-n)/n^ ^ n > 3 , 

+ *» = 2 - (5 ' l} 

We say that U(x,t) is a solution of the Zzmzi problem if for every f G [0,T], U(-,t) 
satisfies the following obstacle problem: 

w G K-t, 

a(w,(f>) > {-L,4>), V0GV, (5.2) 
^ a(w,Tpw) = {—L,ipw), Vtp G W, 



where 



<Lg f|fl' 1 (M n \B e )nC(M n \ J B e ) : 

I £>0 



99 > 0, lim 



M-S-0 U AL {x,t) 



V = {(j) G i/ 1 (M n ), > 0, (/> = on B e (0) for some e > 0}. (5.3a) 

and 

W = VnC 1 (R n ). (5.3b) 

Theorem 5.1. For given A > and L > 0, U At L(x,t) defined in (|5.ip is the unique 
solution of the problem (|5.2|) 

Proof. We first verify that £/a,l is a solution. It can be found by integrating the function 
V A ,L(x,t) from Lemma l4T5j 

K ' Jo 1 ' \Sl {x) V{x,t)dt, t>s(x) 
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where 

s (x) = p^(\x\). 

It is clear that AU = L in U > and U = f^j = on T t (U), i.e. when \x\ = p(t). 
A straightforward application of the Green's theorem then yields that U satisfies the 
variational equality and inequality in ()5.2|) . 

To prove uniqueness, we will use comparison for an obstacle problem and for this we 
need the following observation: 

Lemma 5.2. If w satisfies the obstacle problem (|5.2p . it also for every a > satisfies 
the obstacle problem 

an(w, if — w) > (-L, p — w)q for all p £ IC, 

on tt = R n \ B a , where 

JC = {p G H 1 ^) :<p>0, <p = w on dB a ). 

Proof. Fix p G K. Pick < e < a. Find V 6 C^R 71 ), < V < 1 such that V = on 
B a _ £ (0) and ip = 1 on £ a (0) c . Then ^ eW. Define 



V? on5 a (0) c , 
w otherwise. 



Since p\ 9Ba (o) = MdB a (o), P G # \ B a -e). Define 

= (if) - l)w + <p G if^JR"). 

Clearly > and = for \x\ < a — e. Therefore <p G V . We have from problem (15.21) 

aa(w, p — w) = a(w, <p — ipw) > (-L, 4> — ipw) = {—L, p — w) n . 

This holds for every ip G /C. □ 

The main tool is the following lemma: 

Lemma 5.3 (Comparison for the limit problem). Let w\, W2 be two solutions of the 
obstacle problem (15. 2p for some t > with A\, L\ resp. A2, L 2 . If < A\ < A% and 
Li > L 2 > then 

Wi < 1^2 for all x 7^ 0. 



Proof. Let e 



> 0. If n > 3, there exists a > such that 

u> 2 (2;) 



3 



1 2-n 



Ait 



< e and 



|2-n 



A 2 t 



< e for all \x\ < a. 



We can replace \x\ n by — | log |a;| in dimension n = 2. In particular, ioi(:r) < it^O* 1 ) m 
\x\ < a. Lemma l5.2l implies that w\ and W2 satisfy the obstacle problem on S] = B a (0) c , 

a(wi, p - Wi)n > (-Li, p - Wi) n Vp G /Q, i= 1, 2, 

where /Q = G H l (£i) : p > 0, 93 = it;, on |x| = a}. Now we can use the comparison 
for the obstacle problem, see Corollary 5.2, chapter 4 in [15], and that gives 



in 0. 



□ 
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Now we can finish proof of Theorem 15.11 

Lemma 15.31 implies that for any e > we can compare 

U A - £ ,l(x, t) <U (x, t) < U a+£ ,l(x, t), , 

since UA,L(x,t) are solutions of the limit problem (|5.2p . 
After taking limit e — > we conclude that 

U{x,t) = U A , L {x,t). 

□ 



6 Uniform convergence of u x 

The following lemma was proven in [8]: 

Lemma 6.1 (cf. pH Lemma 4.1]). For given g satisfying (jl.ip and (jl.2p . there exists 

a constant, denoted such that if C M n is a bounded measurable set and if 

{u e } £> C L 2 (Q) is a collection of functions such that u £ — > « strongly in L 2 (f2) as 
e — >• 0, then 

lim / — — ; -u 6 (x) dx = i ( -)u(x) dx a.e. uj (6.1) 

s-+oJ u g(x/e,u) Jq\9/ 



Recall that in section 2, we fixed u for which (|6.ip holds. The goal of this section is 
proving the following theorem: 

Theorem 6.2. The functions u x converge to Ua,l as A — ?■ oo locally uniformly on R n \ 
{0} x [0, oo), where Ua,l is the unique solution of the limit problem from Theorem 15 .1\ 

with A = C* from Lemma 14.41 and L = /|\ from Lemma 16. II 

Proof. We need to prove uniform convergence on sets {\x\ > e, < i < T}. Hence fix 
T > and e > 0. The functions w A (-, t) satisfy the obstacle problem (|3.4p for each t > 
and A > 0. 

Using Lemma H~2"l we can find P2 > and to > such that 

fi t (u) C 5 p2 ^ w for all t > t (6.2) 

and 

flt(u) C £ P2 for all < t < t 
Inclusion (16. 2D is preserved under rescaling when n > 3. In dimension n = 2, we have 

4 (u A ) C B iz(xt) for all t > t /X. 

p2 -R(A) 

But 

ft(At) ^ 7l(AT) /TlogA\ 1/2 mm 

— -. - ' < ~~ttt~ ~ t^t -> r i/2 as A -> oo. 

ft(A) - 7£(A) V lo g AT / 
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Thus we can find p > such that tt t (u x ) c!l:= B^(0) for all < t < T and A > 1. 
Define Q £ := {x G fi : |x| > e} and Q e := ft e x [0, T], 

Now find Ao > 1 such that K x C B £ / 2 (0) for all A > Ao- Since u x (-,t) is a solution of 
the rescaled obstacle problem (|3.4p . we can use the regularity estimates for the obstacle 
problem from [15} Theorem 5.2.4] to conclude that 



Au x (-,t) 



Li(Sl e/2 ) 



< 



1 



for all 1 < q < oo. 



Let p+ be the radially symmetric supersolution from the proof of Lemma l4.2i Since 
the rescaled p\ converges uniformly on the set Q £ / 2 to the function Va,l as A — )• oo for 
some A and L by Lemma 14.51 we can find a constant C\ such that p\_ < C\ on Q £ /2 for 



all A > Aq. Now recall that due to Theorem 



we can express u x as 



u x (x,t) 



[ v x (x,s)ds< [ pl(x,s) ds <C X T for (x,t) G Q £ / 2 - (6.3) 
Jo Jo 



In particular, 



is bounded uniformly in t G [0,T] and A > Aq- 



Thus we can use the standard elliptic regularity results (see [9], for instance) to find 
constants < a < 1 and C2, independent of t G [0, T] and A > Aq, such that 



U X (;t) 
U X (;t) 



<C 2 , 



for all < t < T, A > A . 

<C 2 . 

C*°. a (Q £ ) 

Using (|6.3p again, we have |ii A (x, t) — u x (x, s)\ < C3 |i — s\ and we conclude that 

,A 



u 



< Ci{C 2 , C3), 



for all A > Ao- 

Using the standard diagonalization argument and Arzela-Ascoli, we can find a sub- 
sequence u Xk that converges locally uniformly on sets Q £ , e > 0, to a function u as 
k — > 00. Due to the compact embedding of H 1 in H 2 , the uniqueness of the limit and 
the bound on the i? 2 -norm implies that also u Xk (-,t) — > u(-,t) in i? 1 -norm on sets Q £ 
as A — > 00 for every < t < T, e > 0. 

Finally, in the two following lemmas, we will show that u is a solution of the limit 
problem (|5.2|) . Using the uniqueness of the limit problem, Theorem I5.1} we conclude 
that the convergence is not restricted to a subsequence and we have 



u 



A 



00, 



locally uniformly on W 1 \ {0} x [0,oo), which concludes the proof of Theorem 
Lemma 6.3. W = u(-,t) satisfies 



a(w,<f>) > (-L,4>), 
a(w,tpu) = (—L,ijjw), 



V<p G V, 
VV> G W, 



for each < t < T, where L 



defined in Lemma 16.11 and V, W were defined in 
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Proof. Fix t G [0, T] and </> G V. We will denote w k = u Xk (x,t). Then there is e > 
such that (j) = in B £ (0). There is also k$ such that if k > fco, we have n o fc C -B £ (0). 
Define sequence ip k = (j) + w k . Clearly <fi = on Tq* 1 if > ko and thus G K-f k . Since 
u Xk satisfies (|3.4f) . we have 



a(u/\ 0) = a(w k , <p k - w k ) > 



The mapping w t— > a(w, <ft) is a linear functional on H 1 and hence we have a(w 



a(w, 4>) as k — > oo because w 
(—L,(j)). Hence we conclude 



k _v 7777 ffl 



-fT 1 -strongly. Lemma 16.11 also yields 



a(w, <f>) > (-L, 



Now fix tp G VF such that < ^ < 1- There is B e (0) in which ijj = and &o as above. 
Define = (1 — ip)w k . Again ip^ G /C^ fc if fc > ko- Now 

a(w k ,ipw k ) = -a(w k ,(f k - w k ) < -/ — ^>V>* ~~ 10 / = ^ — aP^^V 



The fact that w k — > to L 2 -strongly implies — > ^ii; in L 2 and by Lemma 16.11 

1 



.9 



, V>w ) -> (—L,ipW). 



Due to the lower semi-continuity of the map u> i— >■ a(w,tpw) in i? 1 , we also obtain 

a(w,ipw) < lim inf a(w k ,ipw k ) < lim ( T—,ipw k \ = (—L,i/m). 

fc— s-oo fc— »oo \ g ,A * / 

Finally, since ipw S 7, we have 

a(w,ipw) < (—L,ipw) < a(w,ipw). □ 

Lemma 6.4. We have 



lim 



u(x, t) 



M->o Uc„L(x,t) 
for every i > 0, where C* is the constant from Lemma 14.41 

Proof. We will use comparison with the sub- and supersolutions from Lemma 14.11 

Let C* be the constant from Lemma [4.41 Fix e > 0. Then by Lemma [4.41 there exists 
a > big enough such that 



P{x) 



,2-n 



< 



for larl = a. 



The set {\x\ = a} is a compact subset of W 1 \ K and by the near field limit, Theorem 
.31 there is to > large enough so that 



v(x,t) P(x) 



Therefore we have 



,2-n 



v(x, t) 
n 2—n 



,2-n 



c* 



<2 



for \x\ = a,t > to. 



< e for \x\ = a,t > to. 
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Let p + and p- be the radially symmetric supersolution, resp. subsolution, satisfying the 
boundary conditions 
P± 



a 2-n 



C* ± e on \x\ = a, 



b + = max 6_ = min \x\, 

Then by comparison 

p~(x, t — to) < v(x, t) < p+(x, t — to) for \x\ > a,t> to. 
Lemma 14.51 yields 

p± -)■ V± := Va t ± e ,L ± as A oo, 

locally uniformly on {\x\ > 0, t > 0}. 

Using the formula (|3.5[) for u x we have for A 1//n |x| > \a\ and A > to/t 

u x (x,t) = [ v x {x,s) ds 
Jo 

"t /X 



< j v x (x,s)ds + I p+(x,s - t /X) ds 

JO Jto/X 
to/A rt-to/X 

v x (x,s)ds + / p x _(x,s)ds 



The first term — > as A — > oo since v x (x, s) < C\x\ 2 ~ n by Lemma 14.21 The locally uni- 
form convergence of p+ then implies that the second integral converges to J* * V+(x, s) ds 
as A — > oo by the bounded convergence theorem. The same argument can be used to 
bound u x from below using p_ and we finally have 



/ V-(x, s) ds < liminf u x (x,t) 
Jo 



< u(x, t) < limsupu A (x, t) < / V + (x,s) ds. 

A— ¥oo JO 



(6.4) 



for all (x,t) G (R n \ {0}) x [0,oo). 

Now suppose that n > 3. The integrals on both sides can be found explicitly as in 
the proof of Theorem 15.11 and 

/ V±(x,s) ds = U±(x,t) := U Ct ±e,L ± (x,t), 
Jo 

with 

U C ,±e,L ± {x,t) 

hm — -— = (C* ± e)t. 

|ar|->0 \x\ 2 ~ n 

Then we divide (j6.4j) by |x| 2_n and take the limit |x| — > 0, which leads to 

(C* - e)t < liminf j—jir^— < limsup < (C, + e)t. 

\x\^o \x\ 2 - n ^i^o \A 

We obtain an analogous result in dimension n = 2, after replacing |x| 2_n by — log 

The proof is complete because e > was arbitrary. □ 

This finishes the proof of Theorem 16.21 □ 



N. Pozar 



7 Convergence of v and the free boundary 

In this final section, we want to show the locally uniform convergence of v A and r(t» A ) 
as A — > oo. Define the half-relaxed limits 

v*(x,t) = limsup v x (y,s), 

(y,s),X— >(x,t),oo 

v*(x,t) = liminf v x (y,s), 

(y,s),\—t(x,t),oo 

in {\x\ ^ 0, t > 0}. Let V(x,t) = Vc*,L{x,t) to be the radially symmetric solution with 
delta source at from Lemma 14.51 where C* is the constant from Lemma 14.41 and 
L = (l). Clearly 

n t (V) = {x : < \x\ < p(t)}. 

And by continuity, V* = V = V* . 

Our goal is proving a theorem similar to JSJ Theorem 4.4]. 

Theorem 7.1. {r(w A )} A converges to T(V) locally uniformly with respect to the Haus- 
dorff distance and v x converges locally uniformly to V in (W 1 \ {0}) x [0, oo) as A — > oo, 
and furthermore 

v* = v* = V. 

We will use the ideas developed in [8] . The proof is somewhat less technical because 
V, the limit of u A , is smooth. 

First, we will collect some necessary technical results in the spirit of [8]: 

Remark 7.2. (a) Since v is a viscosity supersolution, v x (-,t) is superharmonic in 
r2i(f A ) and therefore also v*(-,t) is superharmonic in Q t (f*)- 

(b) u(-,t) is a weak solution of — Aw < for every t > and thus v(x,t) > ju(x,t). 

(c) The construction of the viscosity solution v implies that v x is subharmonic in (K^ c 
(see [U Proof of Theorem 3.1, step 2]). Therefore v* is subharmonic in R n \ {0}. 

Lemma 7.3 (cf. [SJ Lemma 4.6]). Suppose (x^,^) G {u Xk = 0} and (xk,tk,Xk) 
(xoo, too, oo) with xq 7^ 0. Then the following holds: 

(a) (xq, to) G {U = 0}. 

(b) If x k e r tk (u x *) then xoo G r tao (u). 

Proof. See Lemma 4.6 in 0. □ 

We need to identify the singularity of v* and v* at x = 0. 

Lemma 7.4. v* and v* defined above have a singularity at with 

,. v*(x,t) v*(x,t) 

lim — ; = 1, hm — ; = 1 

\x\-+o+ V(x, t) |x|->o+ V(x,t) 

for each t > 0. 
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Proof. Fix e > 0. Using comparison with radially symmetric subsolution p_ and super- 
solution p + as in Lemma 16, 41 with 

P±(x,t) n , , | | 

— t, = ± £ lor \x\ = a, 

a z ~ n 

we deduce that there are a > and to > such 

p~(x, t — to) < v(x, t) < p+(x, t — to), for all t > to, \x\ > a 

Since 

p±(x,t - to/A) ->■ Pb*± e ,L ± as A — » oo, 
locally uniformly on the set (R n \ {0}) x [0, oo). The uniform convergence yields 

Vc*-s,L-(x,t) < V*(x,t) < V*(x,t) < Vc*+e,L + (x,t). 

We have the conclusion because e > was arbitrary. □ 

Lemma 7.5 (cf. [S\, Lemma 4.12]). There exists a constant C\ = C\{n,M) such that if 
{xo,to) 6 f2(t> A ) and B r (xo) n S7 A = 0, we have 

\ . , Cir 2 
sup V (-,to) > — — • 

BrOro) *0 

Proof. Recall that v x (-,to) > u A (-,io)Ao (Remark IT.2[) and use Lemma 3.3 in [8]. □ 
Lemma 7.6 (cf. [8, Lemma 4.10]). The following inclusion holds: 

n(v) c «(«,). 

Moreover, 

> V. 

Proof. Recall that v(x,t) > ju(x,t) for all i > 0, x 7^ and that the inequality 
is preserved under rescaling (|3.ip and (j3.5|) . As u A converges to C7 uniformly on sets 
{|x| > e} x [0,T], we have 

v*(x,t) > -U(x,t) 

and we see that £l(V) = £l(U) C Cl(v*). 

Recall that v*(-,t) is superharmonic in Oj(u*) for each t (Remark l7.2l) and behaves at 
zero as ~ V(-,t) by Lemma 17741 For fixed e > 0, comparison of ?/*(•, i) and Vc r «-e,L(-, i) 
yields 

v*(x,t) > Vc t -e,L(x,t), for every t > 0, 
because £1{Vc 9 - s ,l) C fi(V). We conclude by taking the limit e — )• 0. □ 

Lemma 7.7 (cf. [8] Lemma 4.13 (ii)]). The following inclusion of boundaries holds: 

T(v*) C T(V). 

Proof. Proof as in [8], with the use of Lemma 17.31 above. □ 
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Proof of Theorem \7.1\ Since 0,f(v*) is bounded for every time t > by Lemma \4.2\ and 
£l(V) is a simply connected set, Lemma 17, 71 also implies that in fact C Q(V). In 

particular, we see that 

Wv*) C WV) C n(Vc.+e,L) for all e > 0. 

Recall that t;*(-,i) is subharmonic in i?" \ {0} for every t > (Remark I7.2p . Since 
lim^i^o = 1 for all t > 0, we see that v*(x,t) < Vc,+ £ ,L(x,t) for every e > 0. 

After sending e — > 0+, we recover 

V(x,t) < v*(x,t) < v*(x,t) < V(x,t) 

and 

r(«,) = r(«*) = r(n 

To prove the local uniform convergence of boundaries with respect to the Hausdorff 
distance, fix < t\ < t% and let us denote 

r A = r(v x ) n{h<t< t 2 }, r°° = r(V) n {h < t < t 2 }. 

We define the ^-neighborhood of a set A C R n x R, 

Eft (A) := {(x,t) : dist((x,i),A) < <5}. 
We need to show that for every 5 > there exists Ao > such that 

T A C c/ 5 (r°°) VA > A . (7.1) 

and 

r°° c Us(T x ) VA > A . (7.2) 

First, suppose that for some 8 > 0, the inclusion (|7.ip fails infinitely often. In other 
words, suppose that there is 5 > and a sequence {A/ C }^1 1 with A& — >■ oo such that 

r Afc n (c/ 5 (r°°)) c ^ 0. 

Accordingly, we choose a sequence of points (xf,,tk) G r Afc whose distance from r°° is 
bounded by 5 from below, 

dist(r°°,(x fc ,t fc )) >5. (7.3) 

Since the union |J A r A is bounded (recall Lemma [4 . 2 f) . there is a converging subsequence 
(xkjjtkj) ( x o,to) as j — > oo. Moreover, xo 7^ by Corollary 14.21 and also t\ < to < *2- 
Hence Lemma ITTBT b) implies that (xo,to) G r([/) = r(V) and therefore (xo,to) e 
But that is a contradiction with (|7.3p . 

To prove the second inclusion, (|7.2p . we start by proving a pointwise result. Suppose 
that there is 5 > 0, a point (xo,io) e an ^ a sequence {Afc}^! with A& — >• 00 such 
that dist((xo, to)> r Afc ) > 8/2 for every /c. That means that there is r > such that the 
spacetime cylinder D r (xo,to) := B r (xo) x [to — r, to + r] is either 

^r(x ,to) C {f Afc > 0}, 

or 

D r (x ,t ) C |> fc =0}. 

At least one of the inclusions holds for infinitely many k. Therefore we arrive at two 
possibilities: 
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(a) In the case of the first inclusion, Lemma 17.51 and Harnack's inequality in _B r / 2 (xo) 
for all t G [to — r, to + r ] yield 

Q r 2 

< sup v Xk (-,t)<C 2 inf v Xk (-,t), 

4l B r/2 (x ) B r/2 {x ) 

where C 2 does not depend on A& and t. We conclude that w* > in £L/ 2 (xo) x 
[to — r, to + r] and that is a contradiction with (xo, to) £ r°° C r(«*). 

(b) In the second case, we clearly have that = in D r (xo,to) but that is a contra- 
diction with (x ,t ) g r°° C r(«*). 

We conclude that for each 5 > and (xo, to) £ r°° there is Aq, depending on (xo, to), such 
that dist((x ,t ),r A ) < 5/2 for all A > A . But if (x,t) G r°° such that |(x,t) - (x ,t )| < 
(5/2, we have 

dist((x,t),T A ) < dist((x ,t ),r A ) + |(x,t) - (x ,to)| < ^ VA > A . 
The set r°° is compact and therefore can be covered by finitely many open sets, 

r°°C [J {(x,t):\(x,t)-(xlti)\<5/2}. 
i<i<io 

We can set Ao = maxi<j<j Aq and we conclude that 

r°° C U S (T X ) VA > A . 

□ 
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